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ABSTRACT 

The two Large Magellanic Cloud star clusters NGC 1805 and NGC 1818 are approximately the 
same chronological age (~30 Myr), but show different radial trends in binary frequency. The F-type 
stars (1.3 - 2.2 M 0 ) in NGC 1818 have a binary frequency that decreases towards the core, while the 
binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with 
a peak in the core). We show here, through detailed Wbody modeling, that both clusters could have 
formed with the same primordial binary frequency and with binary orbital elements and masses drawn 
from the same distributions (defined from observations of open clusters and the field of our Galaxy). 
The observed radial trends in binary frequency for both clusters are best matched with models that 
have initial substructure. Furthermore, both clusters may be evolving along a very similar dynamical 
sequence, with the key difference that NGC 1805 is dynamically older than NGC 1818. The F-type 
binaries in NGC 1818 still show evidence of an initial period of rapid dynamical disruptions (which 
occur preferentially in the core), while NGC 1805 has already begun to recover a higher core binary 
frequency, owing to mass segregation (which will eventually produce a distribution in binary frequency 
that rises only towards the core, as is observed in old Milky Way star clusters). This recovery rate 
increases for higher-mass binaries, and therefore even at one age in one cluster, we predict a similar 
dynamical sequence in the radial distribution of the binary frequency as a function of binary primary 
mass. 

Subject headings: (stars:) binaries: general - galaxies: star clusters: individual (NGC 1805, NGC 1818) 
- (galaxies:) Magellanic Clouds - stars: kinematics and dynamics - methods: nu¬ 
merical 


1. INTRODUCTION 


Binary, and higher-order multiple, stars are ubiquitous, 
and comprise a rel atively large fracti o n of the stars i n star 
forming reg i ons (iGhez et, all 1 19931 iKohler et al.1 120001 : 
iKraus et al.l 120111: iBatel 2012t Sana et all I2013IL open 


clust e rs (iMermilliod et ah 1992t Patience et al.1 1 1 9981 
[2001 IGeller et al.ll2008l 120101: IHole et al.l^OOlT and the 
Galactic field (IRaghavan et all 120101) . For solar-type 
stars, the binary frequencjQ in the Galactic field is ap¬ 
proximately 50%, while for the most massive stars, the 
binary frequency increases to about 70% (and perhaps 
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1 We follow the usual convention and define the binary frequency 
as /b = ATb/(Ar s + Ab + ...), where N^, is the number of binaries, 
N s is the number of single stars, and signifies higher-order 

multiples. 


even higher, e.g. IRaghavan et al.1120101 : ISana et al1l2012t 

iCaballero-Nieves et al.ll20ldTu 

Binaries in the Galactic field live essentially in isola¬ 
tion, and only the very widest field binaries (e.g., those 
with separations over ^10 3 AU) are in danger of encoun¬ 
tering passing stars, or having their orbits chan ged dra¬ 
matic ally by the Galactic tidal field dKaib fe Ravmondl 
I2014T) . However, most st ars (with masses >0.5 M^) form 
in denser environments (ILada fe Lada 20Q3; lEvans et al.1 
2 009 : Bressert et al.l l20~10f )7 Though many young, em¬ 
bedded, star c lusters quickly dissolv e to populate the 
Galactic field (|Adams fe Mverdl2001[f . close stellar en¬ 
counters within these birth environments are capable of 
significantly modifying the orbits of, and even disrupt¬ 
ing, binary systems that would otherwise be stable in 
the Galactic field. Therefore, our interpretation of the 
observed binary populations in star clusters and the field, 
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as well as our understanding of star formation in general, 
relies on how a population of stars and binaries evolves 
through this more dynamically active early stage in a 
star cluster. 

Within star clusters, binarie s are typ ically discussed as 
either being “hard” or “soft” (lHeggielll975D . Hard bina¬ 
ries have high binding energies relative to the kinetic en¬ 
ergies of stars moving throughout the cluster, and there¬ 
fore encounters tend to make hard binaries more tightly 
bound, or harder. The converse is true for soft bina¬ 
ries, and indeed an encounter involving a soft binary can 
completely unbind the system. 

Observations of binary populations in young star clus¬ 
ters are very valuable for our understanding of the early 
phases of star cluster dynamics. Two of the youngest 
well-sampled star clusters that have such observations 
are the Large Magellanic Cloud (LMC) star clusters 
NGC 1805 and N GC 1818. both with an age of approx¬ 
imately 30 Myr (Ide Griis et al.l [2002t ILi et al.l [20 1 .'{ FI. 
For bo th star clusters. lde Griis et al.1 (120131) and lLi et al.1 
(120131) measure the binary frequency as a function of ra¬ 
dius from the cluster center. Interestingly, they find that 
the two clusters show different trends, where the radial 
distribution of the binary frequency in NGC 1818 de¬ 
creases towards the cluster core, and the distribution in 
NGC 1805 is flat or perhaps bimodal (with a peak in the 
core). 

Their result for NGC 1818 is particularly interesting 
because, typically, the binary frequencies of star clusters 
are observed to have the opposite radi al tre nd, in creasin g 
towa r ds the core of the cluster (e.g. iGeller fe Mathieul 
l2012t iMilone et al.1I2012D . This phenomenon is under¬ 
stood to be a result of mass segregation, where, in a 
sample containing binaries with primary stars of similar 
masses to the single stars, the more massive binaries sink 
towards the cluster center, due to two-body relaxation 
and dynamical friction. 

In IGeller et al.1 d2013bl f. we studied NGC 1818 in de¬ 
tail through (V-body modeling. We found that, for a 
cluster born with soft binaries and no radial dependence 
of the binary frequency, the observed trend of the bi¬ 
nary frequency decreasing towards the cluster core can 
be naturally explained through the early disruption of 
wide binaries by close encounters with other stars (on 
approximately a crossing time). This process occurs pref¬ 
erentially towards the denser core of the cluster, which 
explains the lower binary frequency in the core relative 
to the halo. Over a few relaxation times, the binary fre¬ 
quency evolves under the influence of mass-segregation 
effects, to produce a bimodal radial distribution, and 
eventually a distribution that rises only towards the core 
(as is often observed in older star clusters). 

The observed bimodal-like _radial trend in binary fre¬ 
quency for NGC 1805 (jLi et al.l l2013f ) is similar to the 
second phase in this evolutionary sequence, where rnass- 

2 Ide Griis et alj (120021 ) estimate an age for NGC 1805 of log(t 
[y r ]) = 7.00j)g'j l Q and for NGC 1818 of 7.40j)g'fg (from a thorough 
review of the prior literature), while ILi et all 120131) find an age 
for NGC 1805 of log(t [yr]) = 7.65 ± 0.10 and for NGC 1818 of 
7.25 ± 0.10. This slight disagreement in the cluster ages likely is 
a result of improved isochron e s (and better field star decontami¬ 
nation) employed by ILi et al.1 (2013). For our modeling we use a 
single age of 30 Myr for both clusters, which is well within a 30- 
confidence interval from each of these age estimates. 


segregation effects begin to become important for deter¬ 
mining the radial distribution of the binary frequency. 
However, NGC 1805 and NGC 1818 are approximately 
the same chronological age, and furthermore, our models 
of NGC 1818 suggest th at even ta king t he extremes of 
the age rang e from[de Gr iis et al.l ( 2002 1 and references 
therein) and ILi et al.l (120131) . cannot alone explain the 
differences in their binary frequency radial distributions 
(with all else being equal between the clusters). 

The two key differences between these clusters, from 
a dynamical evolution perspective, is that NGC 1805 
is less massi ve and more centra ll y con centrated than 
NGC 1818. iMackev fe Gilmorel ( 20031 ) estimate that 
NGC 1805 and NGC 1818 have masses of log(M/M Q ) = 
3.52 ±0.13 and 4.13±g it, respectively, and core radii of 
1.33 ± 0 .06 pc and 2.4 5 ± 0.0 9 pc, respectively. Conse¬ 
quently, Ide Griis et al.l (120021 ) estimate that the current 
half-mass relaxation time in NGC 1805 is 4 to 5 times 
shorter than in NGC 1818 (Ide Griis et al.ll2002l ). 

Here we use (V-body modeling to investigate whether 
the observations of both NGC 1805 and NGC 1818 can 
be reproduced with similar initial binary populations, 
all un dergo ing the evolutionary sequence discussed in 
IGeller etak (2013b). If so, this may suggest that binaries 
can form with similar properties across diffe rent environ¬ 
m ents in the LMC, as is also discussed in IGeller et all 
(l2013al) in the context of open clusters and the field in 
our Galaxy. 

We briefly discuss the setup of our simulations in Sec¬ 
tion [2] In Section [3] we compare the models directly to 
the observati ons o f NGC 1 805. We then expand upon the 
analysis of lGeller et al.l ; 2013 1 Li and study the evolution 
of the binary frequency as a function of binary primary- 
star mass in Section [I] We offer our explanation of the 
different observed trends in the radial dependence of the 
binary frequency in NGC 1805 and NGC 1818, through 
comparisons to our iV-body models, in Section [5] Fi¬ 
nally, in Sections [6] and [T] we provide a brief discussion 
and conclusions. 


2. TV-BODY SIMULATIONS 

We perform a modest grid of d irect jV-body simula¬ 
tions using the NB0DY6 code ( Aarsethl 2003 ) targeted 
at reproducing the observed surface dens ity profile of 
NGC 1805 from IMackev fc Gilmorel (120031) at an age of 
30 Myr. Our modeling procedure here is nearly identi¬ 
cal to that of IGeller et al.1 (|2013al lbl). and therefore we 
point the reader to those papers for additional details 
not discussed here. 

We stress that the initial conditions for the binary pop¬ 
ulations in these NGC 1805 models are nearly identi¬ 
cal to those of our NGC 1818 models from IGeller et al] 
(l2013bl ). Moreover, we input a 100% initial binary fre¬ 
quency (independent of stellar mass and radial posi¬ 
tion) , with binary orbital parameters distributed accord¬ 
ing to observations of solar-type bina ries in young open 
clusters and the Gal actic field (e.g. IGeller et al.1 l201(j 
iRaghavan et al.ll2010D . Most importantly for our later 
discussion, the initial orbital period distribution is log¬ 
normal (as observed for solar-type binaries, with a mean 
of log(P [days]) = 5.03 and er = 2.28), and extends from 
periods of ab out 0.1 days to 10 10 days. As explained in 
IGeller et al.1 (12013H see particularly Figures 3 and 4), 
and below, many of the very wide (soft) binaries are 
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disrupted quickly by encounters with other stars. The 
initial orbital eccentricities are drawn from a Gaussian 
distribution with a mean of e = 0.38 with a = 0.23, as 
is observed for solar-type binaries in the young (~150 
Myr) open cluster M35 and c onsistent with similar bi¬ 
naries in the Milky W ay field dGeller et alJl201(A l2013at 
iRaghavan et al.l[2010f) . We draw binary mass ratios from 
a uniform distribution, but with limits such that the 
mass-ratio is always less than unity and the mass of 
the secondary star is always greater than 0.1 M@, which 
produces a mass-ratio distribution that is approximately 
of the f orm dN/dq oc q ~ 0,4 , within the observe d mass 
range of lLi et all (1201311 : see lGeller et alJl2013bl for fur¬ 
ther details. The only differences that are imposed on 
the initial binaries in these NGC 1805 models from our 
NGC 1818 models arise from the differences in observed 
cluster masses and structures, as we explain below. 

We choose to investigate two different initial virial ra¬ 
tios of Q = 0.5 (equilibrium) and Q = 0.3 (collaps¬ 
ing), and two different degrees of substructure, using 
fractal distributions of degree^ D = 3 (smooth) and 
D = 2 (clumpy). We follow the same procedure as in 
iGeller et al.1 (l2013bh to generate these initial conditions, 
namely smooth distributions are generated within the 
NB0DY6 code and substruc tured initial conditions are gen¬ 
erated using McLuster (iKiipper et al.ll201~l ). In short, to 
set up the clumpy models, McLuster first defines a frac¬ 
tal distribution of stars within a unit sphere, and then 
“folds” this distribution with a given density profile. For 
simplicity, we choose to begin all mod els with an underly¬ 
ing Plummer density profile (lPlummerl 1191 1): substruc¬ 
ture is imposed on top of the Plummer model for those 
simulations with D = 2. Thus, all models begin centrally 
concentrated, as defined by the Plummer scale radius (or 
equivalently, the virial radius, as discussed below). At an 
age of 30 Myr in the simulations, our initially substruc¬ 
tured models relax to h ave sm ooth dens ity distributions 
(e.g., see Figure 1 of lGeller et al.l[2013bl l. 

Likely there are other potential combinations of ini¬ 
tial Q and D that could reproduce the observations of 
NGC 1805 at an age of 30 Myr. We choose here to in¬ 
vestigate [Q,D\ = [0.5, 3] as a baseline model to compare 
against collapsing ( Q < 0.5) and substructured {D < 3) 
models. In general, other studies (including our mod¬ 
eling of NGC 1818) have simulated star clusters with 
Q ranging from about 0.1 to 0.75, where Q > 0.5 are 
initially expanding clusters. Observations suggest that 
many clusters form with subvirial (Q < 0.5) velocities 

dPeretto et al.ll200fiH Andre et al.ll2f)07t iTobin et al.ll2009l 

iProszkow et all [2009 1. and therefore we choose to focus 

on subvirial models and investigate a moderate value of 
Q = 0.3. Likewise, other studies have investigated a 
range in D 1 from about 1 to 3. Again, here we choose a 
moderate value of D = 2 for our substructured models. 
We do not attempt to identify the most probable specific 
initial values of [Q,D\ through our simulations (and in- 

3 We follow iGoodwin L Whitworth) 1120041 1 and IKiipper et all 
{ 20111 1 to define the degree of substructure with the parameter, D, 
such that the probability that a given particle will have a nearby 
“child” is 2 ^ d ~ 3 ), where D = 3 produces no substructure and 
D < 3 produces fractal density dist ributions. We refer t he reader to 
IGoodwin fe Whitworth! 12004 11 and IKiipper et ahl (1201 11 1 for specific 
details on the algorithm used by McCluster for defining the initial 
fractal density distributions 


deed this is likely not possible to any relevant precision, 
given the observations); instead we investigate for differ¬ 
ences between our baseline (smooth, equilibrium) model 
and modestly subvirial and/or substructured models. 

Compare d to NGC 18 1 8, NG C 1805 is relatively com¬ 
pact. The lElson et, ~aI1 (119871 EFF) profile fit to the 
observations of NGC 1805 bv lMackev fc Gilmore! (120031 1 
gives a = 6.84 ± 0.42 arcsec (1.66 ± 0.10 pc, using the 
canonical LMC distance modulus of 18.5, which equates 
to a scale of 4.116 arcsec pc -1 ) and 7 = 2.81 ± 0.10, 
which can be converted into a King core radius of 
1.33 pc. For reference a similar calculation f or NGC 1818 
yields a core radius of 2.45 pc. iMackev fe Gilmorel 
(12003H found a central surface brightness in NGC 1805 
of log^o = 3.49 ± 0.02 L 0 pc -2 (which implies a cen¬ 
tral mass surface density of 155 ± 7 M^ pc~ 2 , given the 
mass -to-light ratio of 0.05 adopted bv lMackev fc Gilmore! 
12003 1. The implied total mass of NGC 1 805 is then 
log(M/M 0 ) = 3.52 ± 0.13. For reference. Ide Griis et al.1 
([20021) found a mass of NG C 1805 of 2.8lg;g x 10 3 M 0 , 
while l.lohnson et al .1 ( 1200 1 ) found a somewhat higher- 
mass of 6000 M@. 

We aim to create a model of NGC 1805 that reasonably 
matches these parameters at an age 30 Myr, by adjust¬ 
ing the initial virial radius (or equivalently, the Plum¬ 
mer scale radius), and initial total mass. From com- 
paris ons to our previous simulations of ot her star clus¬ 
ters dHurlev et al.l2005HGeller et alJ2013al lHl. we predict 
that only a small amount of mass loss will occur over the 
30 Myr timescale of interest. To begin, we ran a series 
of trial simulations of NGC 1805 over a grid in initial 
number of stars and virial radius, and determ ined that 
starting with 7600 stars, drawn from a iKronpd ( 1200 1.) 
initial-mass function (with masses between 0.1 M 0 and 
50 M 0 ) will evolve to be within the observed mass range 
of NGC 1805 at 30 Myr. This initial number of stars and 
mass function results in an initial total cluster mass of 
~4700 M 0 (already within the range of masses estimated 
from the observations of NGC 1805, though towards the 
high end). The true initial number of stars, and corre¬ 
sponding initial total cluster mass, of NGC 1805 is most 
likely not exactly what we have chosen, In general, be¬ 
ginning with more stars (with all else being equal) would 
result in a longer initial half-mass relaxation time, and 
therefore less rapid dynamical evolution. The opposite 
is true for a cluster with a lower number of initial stars. 
As we discuss below, our simulations do indeed reside 
within the observed mass range of NGC 1805 at an age 
of 30 Myr, which is sufficient for our purposes. 

To determine the initial virial radius for a given model, 
we compared results from our initial set of trial simula¬ 
tions to the observed surface density profile at an age of 
30 Myr. We find initial virial radii between 4 and 6 pc (or 
equivalently Plummer scale radii between 2.36 and 3.53 
pc) fit well to the observations at 30 Myr. Specifically 
we will present models with initial virial ratios, fractal 
dimension, and virial radii, [Q, D, Ry (pc)] = [0.5, 3, 4], 
[0.5, 2, 4], [0.3, 3, 4], [0.3, 2, 6 ]. We do not model the 
embedded phase of the cluster here, and instead begin 
our simulations at t = 0 after gas expulsion and with all 
stars on the zero-age main sequence. We compare the 
simulations at 30 Myr to the observations in Section [3] 

To simplify our comparison between models of 
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Fig. 1.— Projected radial mass surface density profile for one 
NGC 1805 simulation comp ared to the EFF pr o file fi t to the ob¬ 
servations of the cluster by IMackev Sz Gilmore! (120031 ). We show 
the same [Q,D] = [0.5, 2] simulation here as in Figure [2] which 
matches the observed radial distribution of the binary frequency 
particularly closely. (All of our simulations reproduce the observed 
surface density profile.) The binned and hatched area shows the 
region from the simulations within which fall 95% of the 1000 ran¬ 
dom sight line s. The solid gray band sh ows the region encom¬ 
passed by the IMackev &; Gilmord (120031) EFF model, with pa¬ 
rameters log^zo = 3.49 zb 0.02 L ,0 pc -2 (155 zb 7 Mq pc -2 ), 
a = 1.66 ± 0.10 pc, 7 = 2.81 zb 0 . 10 . 

NGC 1805 and NGC 1818, we place these simulations 
of NGC 1805 in the same o rbit around the LMC as we 
used in I Geller et al.l (|2013bl) for the NGC 1818 simula¬ 
tions. Specifically we start the cluster at 3.3 kpc from 
the center of a point mass of 10 10 Mq on a linearized 
circular orbit. NGC 1805 is ob served to be at about 3.86 
- 4.00 degrees (~3.4 - 3.5 kpc IMackev fc Gilmorel[20031) 
from the center of the LMC, but little more is known 
about the cluster’s orbit. This, and the very minor effect 
on the cluster from the LMC tidal field over 30 Myr, sug¬ 
gests that this small difference between the observations 
and our choice of orbit is unimportant for the present 
analysis and understanding. We also note that, in addi¬ 
tion to both NGC 1805 and NGC 1818 being at a similar 
distance from the center of the LMC, they are also on 
the same side of the parent galaxy, and therefore are 
also affected in a similar manner by the tidal field of the 
Small Magellanic Cloud (though this minor effect is not 
modeled here). 

We produce 20 simulations at each combination of 
virial ratio ( Q ) and fractal dimension ( D ), drawing from 
the same initial parameter distributions but varying the 
initial random seed, in order to address the stochasticity 
in TV-body simulations Q In total we present results from 
80 simulations of NGC 1805. 

3. COMPARISON BETWEEN OBSERVATIONS AND 
SIMULATIONS OF NGC 1805 

All of the simulations listed in Section [5] evolve to be 
consistent with the observed total mass, total binary fre- 

4 It is a fairly standard practice to run at least 10 simulations, 
with initial conditions drawn from the same parent distributions, to 
assess the ex pected variation among N -body models of a given st ar 
cluster (e.g., IParker & Goodwinll2012l ; IKouwenhoven et all 1 20141) . 
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Fig. 2.— Binary frequency as a function of radius from the clus¬ 
ter cen ter, comparing one NGC 1805 simulation to the observa¬ 
tions of I Li et ahl (120131) . For both the observations and the simula¬ 
tion, all single stars and binaries with primary-star masses between 
1.3 Mq and 2.2 Mq are included. For binaries, we only count those 
with a mass ratio (q = 1712 /Tni, where mi > m 2 ) >0.55, and each 
binary with q < 0.55 is counted as one single star (due to obser¬ 
vational limitations). The observations are shown in colored stars 
with error bars, and indicate the total binary frequency inside the 
radius where the point is located. In the gray-scale distribution, 
we show the results from the same [Q,D] = [0.5, 2] simulation as 
in Figure IT! We use the same cumulative accounting of the binary 
frequency at the same radial locations as the observations, but for 
the simulation, we plot the distribution of the binary frequency 
at each radius resulting from the 1000 different projected sight 
lines (in bins of fo = 0.005, with darker values indicating higher 
likelihood). This simulation very closely reproduces the observed 
bimodal distribution in binary frequency as a function of radius. 

quency and surface density profile of NGC 1805 at an 
age of 30 Myr. As discussed above, NGC 1805 is ob- 
serve d to h ave a mass between 2000 Mq and 6000 Mq 
(and IMackev fc Gilmorel [2003! derived a mass between 
2450 Mq and 4470 M 0 ). At 30 Myr our models ra nge in 
mass from about 3838 Mq to 4312 Mq. iLi et al.l (|2013t l 
found a binary fraction for NGC 1805 between about 
32% and 40% for binaries with primary masses between 
1.3 Mq and 2.2 Mq, mass ratios >0.55 and within a ra¬ 
dius of 45 arcsec (or about 11 pc at the distance of the 
LMC) from the cluster center. Within the same obser¬ 
vational limits, our simulations have binary frequencies 
that range from about 21% to 36% (and total binary 
frequencies between 73% and 82%), at 30 Myr. 

In Figure |T| we show one example comparison to the 
observed surface density profile of NGC 1805 for a par¬ 
ticular simulation that also closely matches the observed 
radial distribution of the binary frequency (see Figure^. 
Here, as we will also do throughout the paper when an¬ 
alyzing simulations in projection, we take 1000 different 
randomly chosen lines-of-sight through our model and 
combine the results. Strictly, the orbit of the simulated 
cluster within the LMC potential defines a true line-of- 
sight projection for the cluster relative to an observer 
on Earth. However, given the uncertainties in the true 
cluster orbital parameters, our simplistic modeling of the 
LMC potential and random asymmetries within the ini¬ 
tial conditions of substructured simulations, we prefer to 
select many random lines-of-sight and analyze the result¬ 
ing distributions as an ensemble. In general we find that 
the results (e.g., surface density profile, radial distribu¬ 
tion of the binary frequency, etc.) are not particularly 
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TABLE 1 

Summary Table of N-body Simulations 


Q 

D 

Ry 

-^sims 

M c 

o- 0 (obs) 

/b(tot) 

/b(obs) 

P 2 (obs) 



(pc) 


(M 0 ) 

(M 0 pc" 2 * * * ) 



(%) 

0.5 

3 

4 

20 

4092 ± 12 

105 til 

0.7637 ± 0.0023 

0-29 + 8 ;gJ 

27.840 

0.5 

2 

4 

20 

4017 ± 16 

148 t^g 9 

0.7450 ± 0.0022 

0.30 ±°;° 4 

59.315 

0.3 

3 

4 

20 

4098 ± 12 

121 tH 

0.7578 ± 0.0024 

0.28 i °;° 3 

34.030 

0.3 

2 

6 

20 

4099 ± 20 

108 tV 4 0 3 

0.8004 ± 0.0024 

0-31 ±8;81 

49.610 


sensitive to the specific choice of projection. 

We also provide summary information for our simu¬ 
lations in Tabic [1] including the initial virial ratio ( Q ), 
fractal dimension ( D ), and virial radius (i?y), the num¬ 
ber of simulations with these initial parameters (Al s i ms ), 
and the 30 Myr total cluster mass (M c ), central sur¬ 
face density (cro(obs), only including stars with V magni¬ 
tudes <_25 J _ag£roximBtelythe faint limit of observations 
from iMackev fc Gilmore! I2003D . total binary frequency 
(/b(tot)), observable binary frequency (/b(obs), fo r stars 
in the simulations within the observed range of iLi et all 
120131 stated above) and the percentage of simulations 
that match t he observed binary frequency radial distri¬ 
bution of lLi et all (120131 ) to within < 3er (P x 2(obs)). For 
all 30 Myr values of a given [Q, D] pairing, we provide 
the means from the 20 individual simulations. For M c 
and /b(tot), we also provide the standard error on the 
mean, while for <7o(obs) and /b(obs) we provide the up¬ 
per and lower limits within which lie 95% of our random 
sight lines across the 20 simulations. 

In Figure [2] we compare the ob served ra dial distribu¬ 
tion of the binary frequency from ILi et al.1 (120131) to the 
same simulation shown in Figure [T] Note the very close 
agreement in both the binary frequency and radial distri¬ 
bution. Furthermore, this simulation reproduces the bi- 
rnodal structure of the observed radial distribution of the 
binary frequency while drawing initial bina ries from the 
same d istributions as our NGC 1818 model dGeller et al.1 
I2013bl) , which has a distribution of binary frequency at 
the same chronological age that falls towards the core. 
We discuss the reasons for this difference in Section [5l 

Looking now at the full grid of NGC 1805 simulations, 
we compare the observed and simulated radial distribu¬ 
tions of the binary freque ncy thr ough a y * 2 analysis, using 
the same procedure as in IGeller et all (12013bl) . For clar¬ 
ity, we repeat the definition here: 

2 \ - (/e(< r ) - / 0 (< r)) 2 

X 4- e[f 0 « r)] 2 5 1 J 

where / e (< r) and / 0 (< r) are the simulated and ob¬ 
served cumulative binary frequency inside the given ra¬ 

dius, respectively, and e[/ 0 (< r)] is the error on the 
observed cumulative binary frequency at the same ra¬ 

dius. We have not constrained the simulation results, 

and therefore we have eight degrees of freedom (one for 
each bin in Figured]). 

The distributions of these y 2 values, combining results 

from each of the 1000 projections for each of the 20 mod¬ 
els with a given [Q,D] pairing, are shown in Figured] Ad¬ 
ditionally, we provide the percent of sight lines through 
these simulations that match the observed radial distri- 



Fig. 3.— Distribution of y 2 values comparing the radial distri¬ 
bution of the binary f r equen cy from the NGC 1805 7V-body simu¬ 
lations to the lLi et al.1 12013 ') observations (shown in Figurc[2]). We 
define the y 2 value as in Equation [T] and perform this analysis for 
all 1000 sight lines to each of our simulations with the given virial 
ratio ( Q ) and fractal dimension (D). The top and middle pan¬ 
els each combine the sets of simulations that either have the same 
fractal dimension (top) or virial ratio (middle), while the bottom 
panel shows each set of simulations separately. The gray vertical 
line marks the y 2 value that indicates a 3cr difference between the 
observations and simulations (at the given number of degrees of 
freedom). 

bution of the binary frequency to within < 3<t in Table [T] 
First, we find that the initially substructured Simula- 
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tions (D = 2) have a distribution of y 2 values shifted 
significantly to lower values than that of the initially 
smooth simulations (D = 3; top panel of Figure [3]). A 
Kolmogorov-Smirnov test comparing these two distribu¬ 
tions returns a probability < 10”' that the two distribu¬ 
tions are drawn from the same parent population. Con¬ 
versely, the distribution of y 2 values for the NGC 1805 
simulations initially in equilibrium (Q = 0.5) cannot be 
distinguished from that of the initially collapsing simula¬ 
tions (Q = 0.3), as shown in the middle panel of Figure[3l 
Thus we find a preference for substructured initial con¬ 
ditions in our models of NGC 1 805 (whic h we als o found 
for our NGC 1818 models in lGeller et al.ll2013bD . 

Interestingly though, only about half of the simulations 
reasonably reproduce the observed radial distribution of 
the binary frequency. The gray vertical line shows the y 2 
value for a 3<r deviation of the simulation from the ob¬ 
servations (at y 2 ed ~ 2.95). We can consider the ~50% 
of simulations with higher y 2 values than this limit to 
be poor fits to the observations. In general the models 
that fail to reproduce the observations do so either be¬ 
cause they have a different total binary frequency or a 
radial distribution of the binary frequency that only de¬ 
creases towards th e core. Thi s seco nd point is in line with 
the results of iGeller et al.1 (|2013bl l. where we find that, 
for smooth initial conditions, the simulations require at 
least one to two half-mass relaxation times to produce 
a bimodal distribution, and NGC 1805 has only lived 
through about one half-mass relaxation time, or less (see 
Section [5]). 

This may indicate that there exists a more accurate 
choice of initial [Q,D], or other initial parameters, that 
best reproduces the observed cluster today. However, our 
goal here is not to identify the exact initial parameters of 
the cluster, but to investigate whether both NGC 1805 
and NGC 1818 could reasonably have been born with 
very similar primordial binary populations. We find that 
we can closely reproduce the observations from a subset 
of our simulations, and therefore we conclude that the 
observations of both NGC 1805 and NGC 1818 can be 
reproduced from simulations beginning from very simi¬ 
lar initial conditions (with appropriate modifications to 
the initial cluster mass and virial radius). In the follow¬ 
ing section, we step back and examine the overall evolu¬ 
tionary trends of the binaries in our different NGC 1805 
models. 

4. DEPENDENCE OF THE RADIAL DISTRIBUTION OF 
THE BINARY FREQUENCY ON STELLAR MASS 

We can reproduce the two different trends in the ob¬ 
served radial distribution of the binary frequency in 
NGC 1805 and NGC 1818 through A-body simulations 
that begin with the same initial binary frequency and 
binary parameters drawn from the same initial distribu¬ 
tions, at the same chronological age of 30 Myr. Moreover, 
we do not need to invoke a different primordial binary 
population, or other changes to the primordial popula¬ 
tion between these two clusters (other than a difference 
in initial cluster mass and virial radius), to explain these 
observations. We begin here to investigate the differences 
between these two clusters, starting with an analysis of 
the b inary frequency as a function of stellar mass. In 
iGeller et al.1 (|2013bl) we investigated the effects of binary 
disruption, two-body relaxation and mass segregation as 



Fig. 4.— Distributions of binary frequency as a function of radial 
distance from the cluster center for binaries with different primary- 
star masses. We cover the same radial domain as the observations, 
though with a slightly different bin size. Each point represents an 
analysis of the stars inside the given radius. However here we do 
not limit by mass ratio, and we take the three dimensional radius, 
rather than projected (as we did in Figure [2] to compare with the 
observations). Different plotting styles and colors identify the dif¬ 
ferent mass bins and are indicated in the top panel. Each point 
shows the mean value from each of the 20 simulations for the given 
[Q,D] pairing. Uncertainties show the standard errors on these 
means (with error propagation in the plotted ratios). As an indi¬ 
cator of the range in uncertainties from all mass bins, and for ease 
of viewing, we only show uncertainties for the lowest and highest 
mass bin (which have approximately a factor of 200 difference in 
the total number of binaries). In the top panel, we show the binary 
frequency as a function of radius at 30 Myr (the age of NGC 1805) 
for our [Q,D\ = [0.5, 2] NGC 1805 simulations, the set of simula¬ 
tions that are most consistent with the binary observations (e.g., 
see Figure |3| and Table [TJ. We normalize each bin by the value at 
the start of the simulation to remove any stochastic fluctuations 
resulting from drawing small numbers of stars from broad distri¬ 
butions in each of the individual simulations. In the middle panel, 
we plot the ratio of these values to those of our [0.5, 3] NGC 1805 
simulations, the set of simulations that are the least consistent 
with the binary observations. Finally, in the bottom panel we plot 
the ratio of the values in the top panel to a similar analysis of our 
NGC 1818 simulations with [Q,D] = [0.5, 2] (the same initial [Q,D\ 
as shown for NGC 1805 in the top panel). 
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a function of time. Examining the binary sample by stel¬ 
lar mass allows us to probe these effects at one age in the 
simulation, since higher mass stars are expected to un¬ 
dergo dynamical relaxation processes more quickly than 
lower mass stars. 

We focus on the NGC 1805 simulations with 
[Q,D\ = [0.5, 2], as these models most closely reproduce 
the observed radial distribution of the binary frequency 
for the observed mass range (see Figures [2] and [3] and 
Table [T]). We divide the model into mass bins of 0.6 M 0 , 
approximately the mean stellar mass in the simulation 
at 30 Myr. At masses greater than 10 times this mean 
mass, the sample size gets too small to continue this bin 
size (without very large uncertainties), and we therefore 
group these highest mass stars together into one bin. 
(The maximum mass of a star in our simulations at 30 
Myr is ^9 M 0 .) Due to the relatively small number of 
stars in NGC 1805, it is unfortunately not possible to re¬ 
trieve useful binary frequency radial distributions from 
the observations as a function of mass, due to the rela¬ 
tively large statistical uncertainties, and therefore here 
we only analyze the simulations. 

We plot the radial distribution of the binary frequency 
in different mass bins in Figure [4] Here we do not project 
the model, and instead use the 3D radial position for 
each star. We also take all binaries into consideration, 
rather than limiting to m ass ratios > 0.55 as we do for 
direct comparisons to the iLi et all ( 20131 ) observations. 
Although we do not explicitly impose any differences in 
the binary population as a function of primary mass, we 
choose to normalize the results at 30 Myr to those at the 
start of our simulation to remove any stochastic fluctua¬ 
tions that result from drawing relatively small numbers 
of binaries from broad distributions (most relevant for 
the higher mass bins, which have smaller numbers of bi¬ 
naries) . 

Focusing on the top panel of FigurelH we point out two 
general trends in the binary frequency that clearly vary 
with mass. First, at an age of 30 Myr, the overall binary 
frequency increases towards higher-mass stars. Second, 
the binaries with higher-mass primaries show more dy¬ 
namically evolved radial distributions. Furthermore, a 
X 2 test comparing the radial distribution of the binary 
frequency for the lowest-mass binaries to that of any of 
the other mass bins returns a distinction at very high 
confidence (> 6cr). 

This first phenomenon can be understood by also look¬ 
ing at the distributions of semi-major axes for binaries of 
different masses. In the top panel of Figure[5]we show the 
mean semi-major axis within bins of the primary star’s 
mass for both the [Q,D] = [0.5, 2] an d [0.5, 3] models. 
As discussed in detail in iGeller et all (2013b ) and also 
relevant here, the widest binaries are dynamically dis¬ 
rupted early in the simulations, which accounts for the 
overall drop in binary frequency seen in Figure [I] Fig¬ 
ure [5] shows that the higher-mass binaries survive out 
to larger separations than the lower-mass binaries (al¬ 
though this is more apparent for the initially smooth 
model, and we will discuss this difference in Section mil . 
As the semi-major axes for all binaries were drawn from 
the same (broad log-normal) distribution, this result in¬ 
dicates that a larger frequency of higher-mass binaries 
than lower-mass binaries are able to survive disruptions 
by an age of 30 Myr. 


The second trend that we point out above may be ex¬ 
pected from two-body relaxation timescale arguments. 
Returning to the top panel of Figure [4] we see that, for 
instance, the binaries with primary masses less than the 
mean mass have a binary frequency that decreases to¬ 
wards the core of the cluster. Moving to bins contain¬ 
ing higher and higher mass binaries, the distributions 
increase more and more towards the cluster center, and 
also begin to display a bimoda l distribution . _We dis¬ 
cussed this behavior in detail in IGeller et al.l ( 2013bt f as 
a typical trend that a population of binaries may go 
through as a function of time. Here we note that even 
at one given time, the same trend can be observed by 
examining binaries of different masses. Moreover, this 
shifting from a decreasing to increasing trend of binary 
frequency towards the cluster core depends on the re¬ 
laxation timescale, which decreases towards higher-mass 
stars. 

4.1. Differences in the Evolution of the Smooth and 
Substructured Models 

Our models with substructured initial conditions (D = 
2) reproduce the observations more closely than those 
with smooth initial conditions (see Figure[3]and Table[l|). 
Below, we point out a few key differences between these 
two sets of simulations which may help to explain the rea¬ 
sons behind this preference. For simplicity, we will focus 
on the simulations that begin initially in virial equilib¬ 
rium (Q = 0.5). 

We start with the middle panel of Figure [2 where we 
plot the ratio of radial distributions of the binary fre¬ 
quency from the [Q,D\ = [0.5, 3] (initially smooth) sim¬ 
ulations over that of the [0.5, 2] (initially substructured) 
simulations. For nearly all mass bins, this ratio decreases 
towards the cluster core. Moreover, this ratio decreases 
more strongly towards the core as we consider binaries 
with higher mass primaries. A y 2 test distinguishes the 
distribution for lowest-mass binaries from that of the 
highest-mass binaries at > 4.5er (and at > 3<r for the 
highest-mass bin compared to the 0.6 - 1.2 M 0 bin, the 
1.8 - 2.4 M 0 bin and the 2.4 - 3.0 M 0 bin). In other 
words, at 30 Myr, the binary frequency in the core for 
a given binary primary mass is lower in the simulations 
with smooth initial conditions than for those with sub¬ 
structured initial conditions, and this difference is more 
pronounced for higher-mass binaries. 

It may be tempting to conclude that the lower binary 
frequency at 30 Myr in the core of the initially smooth 
models is due to the more efficient disruption of binaries 
in those models, as compared to the initially substruc¬ 
tured models. However, turning to Figure [5j where we 
compare the mean semi-major axes for binaries in bins of 
increasing primary mass at 30 Myr, we see that in fact 
the opposite is true. Again, we began all of our simu¬ 
lations by drawing binaries from the same initial distri¬ 
bution of semi-major axes (with no dependence on ini¬ 
tial masses) and allow dynamical encounters to naturally 
disrupt wide binaries. The simulations that begin with 
initially substructured density distributions are more ef¬ 
fective at disrupting binaries, of all primary masses. At 
an age of 30 Myr, the mean semi-major axis for binaries 
within one initial half-mass radius (r^(0) ~ 3.07 pc) for 
the 20 simulations initially with smooth density distribu¬ 
tions ([Q,D\= [0.5, 3]) is 190 ± 3 AU, compared to the 
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160 ± 3 AU for the simulations with initially substruc¬ 
tured density distributions {[Q,D\= [0.5, 2]) within the 
same radial domain. A Kolmogorov-Smirnov test com¬ 
paring these two semi-major axis distributions returns a 
probability of 2xlO~ 10 that the two are drawn from the 
same parent distribution. 

This distinction is also seen in the top panel of Fig¬ 
ure [5j where the mean semi-major axis in every mass 
bin for the [0.5, 3] simulation is higher than that of the 
[0.5, 2] simulation, and in the bottom panel, where the 
ratio of these values is always greater than one. The 
bottom panel of Figure [5] also shows that the ratio of the 
mean semi-major axes for these two sets of simulations is 
largest for binaries with the highest-mass primaries. This 
result indicates that in the initially substructured simu¬ 
lations, the higher-mass binaries undergo more energetic 
encounters than in the initially smooth simulations. 

These effects can all be understood under the same 
physical framework, as follows. The initially substruc¬ 
tured simulations undergo more rapid relaxation than 
the initially smooth simulations early in the cluster evo¬ 
lution. This situation p hysically resembles the condi¬ 
tions of violent relaxation flLyndemBe ll 1967; Elso n et al.1 
119871 ISpera fe Capuzzo-Dolcettal 1201 51) . where the re¬ 
sponse of particles to processes like mass segregation is 
stronger if there are significant gradients in the back¬ 
ground potential, as we also have here in the presence 
of substructure. This causes the higher-mass binaries 
to segregate more quickly towards the core, producing 
the effect seen in Figure [U The result that star clus¬ 
ters with initially substructured density distributions 
undergo mass segregation more rapidly is well known 
from IV-bo dy star c l uste r models, and was first discussed 
bv iMcMillan et ah (120071! land soon therea f ter by , e.g., 
lAllison et al.1 12009 and TMoeckel fc Bonnelll 2009 using 
more extensive models), as a potential explanation for 
very young star clusters that are observed to have sig¬ 
nificant mass segregation, which is inconsistent with ex¬ 
pectations from standard two-body relaxation processes 
(from unsegregated initial conditions). We confirm this 
result here, now from the viewpoint of binary stars. 

Furthermore, because of pockets of initially very high 
density in the substructured simulations relative to the 
smooth simulations, the semi-major axis distributions for 
binaries of all masses are shifted to lower values, by a 
factor of about 1.2, on average, inside of Th{ 0) (see also 
the top panel of Figure [5|. Additionally the higher-mass 
binaries in the substructured models segregate towards 
these high density regions most efficiently, and therefore 
show an even larger difference in the semi-major axis 
distributions from the initially smooth simulations than 
do the lower-mass binaries (bottom panel of Figure [5]). 

5. A COMPARISON TO NGC 1818 

As presented in iLi et ah] (I2013H and discussed in Sec¬ 
tion [T| in this paper, unlike NGC 1805, the radial dis¬ 
tribution of the binary frequency in the LMC star clus¬ 
ter NGC 1818 decreases (only) towards the core. This 
is intriguing in its own right because older star clus¬ 
ters are consistently observed to have the opposite radial 
trend in binary frequency, but is even more interesting 
when compared to NGC 1805, which is also located in 
the LMC and has the same chronological age, but in¬ 
stead shows a bimodal distribution in binary frequency 



Fig. 5.— Distributions of mean binary semi-major axis for bins 
of different binary primary mass, comparing our NGC 1805 sim¬ 
ulations with [Q,D] = [0.5, 2] and [0.5, 3] at 30 Myr. We only 
include binaries with main-sequence primary stars within one ini¬ 
tial half-mass radius (~3.07 pc) of the cluster center. The points 
show the mean values from each of the 20 respective simulations 
for the given [Q,D\ pairing, and the vertical error bars show the 
standard errors on the means (with error propagation in the bot¬ 
tom panel). The horizontal error bars indicate the bin size; the last 
bin extends to the most massive stars in the cluster (~9 Mq), but 
is truncated here for ease of viewing. In the top panel we compare 
the distributions for the two sets of simulations directly (with the 
D =2 simulations in black and the D =3 simulations in gray). In 
the bottom panel we take the ratio of the results from these two 
simulations. 

with radius. Timescale estimates from the observations 
indicate that the current half-mass relaxation time in 
NGC 1805 is 4 to 5 times shorter than in NGC 1818 
(jde Griis et al.l[2002l l. We expand on this discussion be¬ 
low, by comparing our simulations and through simple 
analytic timescale estimates, to argue that NGC 1818 is 
less dynamically evolved than NGC 1805, and this can 
explain the difference in the radial distributions of their 
binary frequencies. 

We re-analyze the two \Q ,D\ = [0. 5, 2] simulations of 
NGC 1818 from I Geller et all ( 2013bl) . and compare the 
results to those of our 20 [0.5, 2] NGC 1805 simulations 
in the bottom panel of Figure [1] Here we do not plot 
the distributions of the highest-mass binaries because the 
uncertainties become so large for the NGC 1818 simu¬ 
lations (due to small sample sizes) that a precise com¬ 
parison is not possible. First, we find that the overall 
binary frequency in NGC 1818 is lower. This is due to 
the higher total cluster mass of NGC 1818, and there¬ 
fore the higher velocity dispersion, which results in more 
binaries being disrupted (see also iSollimal 12008 ). More 
importantly for this discussion, the higher-mass binaries 
in NGC 1818 have a much lower binary frequency to¬ 
wards the core than in NGC 1805, while there is no such 
radial trend for the lower-mass binary comparison. For¬ 
mally, the distribution for highest-mass binaries can be 
distinguished from those of both the two lowest-mass bi- 
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nary bins at > 3er confidence. We interpret this in an 
analogous way to our interpretation of the middle panel 
of the same figure; NGC 1818 is less dynamically evolved 
than NGC 1805 at the same chronological age. 

This difference in dynamical age is also apparent 
when we compare the radial mass distributions in the 
NGC 1805 and NGC 1818 models in Figure [G] At the 
same chronological age of 30 Myr (and for the same [Q,D\ 
= [0.5, 2] initial conditions), the radial mass distribution 
rises more steeply towards the core in the NGC 1805 
model than in the NGC 1818 model. Moreover, the NGC 
1805 model achieves a higher degree of mass segregation 
than the NGC 1818 model, within the same cluster life¬ 
time. 

This result agrees with simple timescale arguments. 
For instance, we calculate the initial half-mass relaxation 
time for stars with masses equal to the mean mass of 
an object (single or combined mass of a binary) in our 
NGC 1805 [0.5, 2] model to be about 50 Myr, while for 
the NGC 1818 [0.5, 2] mo dels a si milar calculation yields 
about 380 Myr (following iSoitzerl I1987L though the for¬ 
mula does not account f or substr ucture). This result 
agrees roughly with lde Griis et ah! ( 2002 ). who estimate 
from observations that the present-day half-mass relax¬ 
ation time in NGC 1805 is 4 to 5 times shorter than in 
NGC 1818. Interestingly, a simple estimate of the cross¬ 
ing time for these two models yields about 5 Myr and 
8 Myr, respectively. Thus, the early phase of dynamical 
disruption, which creates a decreasing trend in binary 
frequency towards the core, likely operated on a similar 
timescale in both NGC 1805 and NGC 1818, but the re¬ 
covery of the binary frequency in the core through mass 
segregation processes operates ~ 5 - 10 times faster in 
NGC 1805 than in NGC 1818. 

We also note that although NGC 1805 may be more 
dynamically evolved, NGC 1818 still exhibits a similar 
trend in the radial distribution of the binary frequency 
as we find in our NGC 1805 models when examining bi¬ 
naries of different mass bins. In Figure [7] we show the 
result from re-analyzing one of the iGeller et all (2013b) 
NGC 1818 simulation with [Q,D] = [0.5, 2] (top), com¬ 
pared to the observations of NGC 1818 in the same mass 
bins. (Unlike NGC 1805, in NGC 1818 the larger number 
of stars enables us to divide the sample in different mass 
bins, though the uncertainties are large.) Again, we see 
that the lowest mass binaries show a decreasing binary 
frequency towards the cluster core, while moving towards 
higher-mass binaries shows the trend shifting towards an 
increasing binary frequency towards the core (in both the 
observations and simulations). 

Finally, the result for NGC 1818 shown in Figure [71 
may help to explain the differ ences between the obser- 
vations of iLi et, all (l2013fl and lElson et,~aLl (|1998l) , who 
found opposite t rends in the_binary frequency as a func¬ 
tion of radius. ILi et al.l (120131) included binaries with 
primary masses between about 1.3 M 0 and 2.2 M 0 , and 
found a decre asing binary freque ncy towards the cluster 
center, while lElson et all ( 1998 ) included binaries with 
primary masses between 2 M 0 and 5.5 M 0 , and found an 
increasing binary frequency towards the cluster center. 
We include these mass bins in Figure[7](both for binaries 
with mass ratios > 0.55), and find different r adial trends. 
We reproduce the observations of ILi et all ([20 13;. as in 
IGeller et al.ll2013bh for the lower-mass sample, while the 



Fig. 6. — Mass distributions as functions of radius from the 
respective cluster centers in the NGC 1805 (black circles) and 
NGC 1818 (blue triangles) models. Each point shows the mean 
mass inside of the given radius divided by the mean mass of all 
stars inside of two times the initial half-mass radius, r^(0), at a 
time of 30 Myr in the respective model. For both NGC 1805 and 
NGC 1818, we use the [Q,D] = [0.5, 2] simulations (as we also show 
in other figures in this paper). The uncertainties for the NGC 1805 
points show the standard errors on the means across the 20 indi¬ 
vidual simulations. For the uncertainties on the NGC 1818 points, 
we simply show the range in values from the 2 individual simula¬ 
tions. In constructing this figure, we consider each binary as one 
object with a mass equal to the sum of the individual component 
masses. 

higher-mass binaries of the lElson et al.l (Il998f) sample 
show a flat, or perhaps marginally increasing, trend in 
binary frequency towards the core (simil ar to t he con ¬ 
clusions from the N -bod y models of lElson et all 1199811 . 
Ide Griis et all ([20131 1 and lLi et, al.l (1201 .'1 speculated that 
this difference in binary mass may explain the discrepant 
observations, and we bolster th is argument here. We note 
that, as stated in Geller et al.l (l2013bil . our NGC 1818 
models do n ot p ro duce as high a binary frequency as 
observed bv lElson et all ( 19981 using their observational 
constraints in our analysis of the simulations). Therefore 
we do not perform a more detailed comparison to the ob¬ 
servations. We simply offer this discussi on as a possible 
way to reconci le the analyses of lElson et al.l (|1998T) and 
ILi et all (Uni). 

6 . DISCUSSION 

The two rich LMC star clusters NGC 1805 and 
NGC 1818 are both approximately the same chronolog¬ 
ical age, but are observed to have two different trends 
in the radial distr ibution of their binary frequencies. We 
found in IGeller et al.l ( 2013b l that the radial distribu¬ 
tion of the binary frequency can evolve over time from 
one that decreases towards the core (like in NGC 1818) 
to a bimodal distribution (like in NGC 1805) and even¬ 
tually to a distribution that increases only towards the 
core (as observed in many older star clusters). The pro¬ 
cess occurs on a two-body relaxation timescale, and be¬ 
cause NGC 1805 is less massive and more compact than 
NGC 1818, the ex pected rel axati on t ime is shorter in 
NGC 1805. Indeed Ide Griis et all (|2002l) estimated from 
observations that the present-day half-mass relaxation 
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Fig. 7. — Comparison of the radial distributions of the binary 
frequency in different mass bi ns fo r on e si mulation of NGC 1818 
with [Q,D] = [0.5, 2] f rom I Geller et all I2013bl . top) and the obser¬ 
vations from I Li et al.l (I2013I . bottom). The different plotting styles 
for the different mass bins are indicated on each plot. As in Fig¬ 
ure^ we show the binary frequency inside the given radius. We 
only count binaries with mass ratios q = rri 2 /m\ > 0.55, and any 
binary with q < 0.55 is counted as a single star. Fo r th e simula¬ 
tion, we use the same 1000 sight lines as in IGeller et al.l ( 2013b '); 
the points in the top panel show the means of these values, with 
uncertainties (visible when larger than the plotted symbols) in¬ 
dicating the standard error on these means. Uncertainti es on the 
obser vations (derived using Poisson counting statistics, see lLi et al.l 
120131 ) are much larger, and are shown in the bottom panel for two 
representative samples, for clarity. 

time in NGC 1805 may be 4 to 5 times shorter than 
in NGC 1818. We therefore conjectured in iGeller et al.1 
(I2013bl ) that if the binaries in both clusters follow the 
same evolutionary sequence, the different observed ra¬ 
dial trends in the binary frequencies between these two 
clusters may be further evidence that NGC 1805 is more 
dynamically evolved than NGC 1818. 

We confirm this suggestion here through detailed N- 
body modeling of NGC 1805. We can reproduce the ob¬ 
servations of NGC 1805 through (V-body models that 
draw from the same initial binary distributions, frac¬ 
tal dimension, and virial ratio as our NGC 1818 models 
(which reproduce the observations of that cluster), with 
the only differences being the initial masses of the clus¬ 
ters and the initial virial radii. Moreover, our simulations 


show that both clusters could have formed with binaries 
drawn from the same parent population, but are observed 
to be different today because of the clusters’ different 
dynamical ages (where dynamical age is the number of 
relaxation times that the cluster has lived through). 

Furthermore, we show above that binaries with differ¬ 
ent mass primaries undergo this evolutionary sequence of 
the radial distribution of the binary frequency at different 
rates, where the higher-mass binaries enter the binrodal 
and mass segregated phases earlier than the lower-mass 
binaries (e.g., top panel of Figured]). This is expected 
from relaxation time arguments, and we show that this 
difference can also be used to compare the dynamical 
ages of clusters, for instance as we do in the bottom panel 
of Figured] This also highlights the importance of exam¬ 
ining binaries in the same mass range when comparing 
observations across clusters. 

As with our NGC 1818 models, we find that our 
NGC 1805 models that initially have some degree of 
substructure (here a fractal distribution with dimension 
D = 2) more often match the observed radial trend in 
binary frequency. We do not attempt to quantify the 
exact degree or nature of this substructure here. Indeed 
it is doubtful that such an analysis would be possible, 
given the available observations, and that the substruc¬ 
ture in our models is erased before an age of 30 Myr. 
Nonetheless, the preference for substructure is encourag¬ 
ing, as this is consisten t with observat ions of star form¬ 
ing regio ns (]Larson||1995HCartwright fe Whitworth| |2004t 
IKraus &; Hillenbrandll2008HSanchez fe Alfaroll2009| b 

Finally, we note again that, although we focus 
on an age of 30 Myr for both the NGC 1805 and 
NGC 1818 models, b oth Ide Griis et all (I2002L and ref¬ 
erences therein) and ILi et al.1 (|2013l) find the clusters 
to have marginally different ages, though these two pa¬ 
pers e ach find a different cluster as being older. [Li et al.1 
(120131 ) . who used improved isochrones and performed a 
more careful field star decontamination, find a nominal 
age for NGC 1805 of 45 Myr and for NGC 1818 of 18 Myr. 
This slight difference in age is in the right direction to 
help increase the expected difference in the radial distri¬ 
butions of the binary frequency in both clusters (though 
if all else were equal between these clusters, this differ¬ 
ence in age on its own would not be enough to produce 
the different radial trends in binary frequency). 

7. CONCLUSIONS 

The radial distribution of the binary frequency in a 
star cluster evolves with time due to dynamical disrup¬ 
tions from close encounters with other stars and mass- 
segregation processes, and can be used to track the dy¬ 
namical age of a cluster. For a star cluster that is born 
with wide (soft) binaries, the early evolution of the bi¬ 
naries is dominated by disruptions, which decrease the 
overall binary frequency, and establish a decreasing trend 
in binary frequency towards the cluster core, on approx¬ 
imately a crossing time. The rich LMC star cluster 
NGC 1818 is observed in this phase of evolution. On 
a two-body relaxation timescale, dynamical friction and 
mass segregation effects take over, causing the more mas¬ 
sive binaries to sink towards the cluster core, which pro¬ 
duces a bimodal radial distribution of the binary fre¬ 
quency. NGC 1805 is observed in this phase of evolu¬ 
tion. Later on, as binaries towards the halo also begin to 
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experience the effects of dynamical friction, the binary 
frequency is transformed into one that only increases to¬ 
wa rds the core, as is ob served in most older star clusters. 
In iGeller et al'l rt2013bl ) we showed that this evolution¬ 
ary sequence can be tracked by looking at binaries of 
the same mass over time, and that the relevant timescale 
is not necessarily the chronological age, but instead the 
number of relaxation times the cluster has lived through, 
which we refer to as the cluster’s dynamical age. 

We show here that the same evolutionary sequence in 
the radial dependence of the binary frequency, can be 
observed at one chronological age for binaries of differ¬ 
ent primary masses. The higher-mass binaries undergo 
dynamical friction and mass segregation processes at a 
faster rate than the lower-mass binaries, and therefore, 
although all binaries are subject to dynamical disrup¬ 
tions early on, the higher-mass binaries begin increasing 
their core binary frequency more quickly. 

Our detailed iV-body simulations confirm that 
NGC 1805 is dynamically older than NGC 1818. Im¬ 
portantly, we can reproduce the observations of both 
clusters by drawing their stellar populations from the 
same parent population (though starting from a differ¬ 
ent total cluster mass and concentration). Furthermore, 
we show that today’s observations of both LMC clus¬ 
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ters can be reproduced by drawing their initial binary 
populations from distributions that are also consistent 
with observed solar-ty pe binaries in the Milky Way field 
(Rag havan et al.ll2010D and o bservations of youn g Milky 
Way open clusters (e.g. M85. IGeller et, al.lI201flf ) . These 
results are consistent with a hypothesis that the bina¬ 
ries in NGC 1805 and NGC 1818 were born with similar 
properties as those in Milky Way clusters, and suggest 
that binaries may form with similar distributions of or¬ 
bital parameters and masses within a variety of different 
environments. We suggest that the radial distributions 
of the binary frequencies in NGC 1805 and NGC 1818 are 
different today simply because we are catching them at 
slightly different stages along a very similar evolutionary 
sequence. 
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